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Abstract 

We present the resuhs of a complete two-loop calculation at order 0{q^) of the nucleon mass in 
manifestly Lorentz-invariant chiral perturbation theory. The renormalization is performed using 
the reformulated infrared renormalization, which allows for the treatment of two-loop integrals 
while preserving all relevant symmetries, in particular chiral symmetry. 
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I. INTRODUCTION 



Chiral perturbation theory (ChPT) [l|, 0] is the effective field theory for the strong in- 
teractions at low energies and has been highly successful in the application to purely mesonic 
processes (for a review see, e.g., Ref. 0]). The extension to the one-nucleon sector was first 
addressed in Ref. and originally seemed to be problematic due to the presence of the nu- 
cleon mass as an additional mass scale that does not vanish in the chiral limit. When using 
dimensional regularization in combination with the modified minimal subtraction scheme 
there is no direct correspondence between the loop expansion and the chiral expansion. The 
first solution to this power counting problem was given by heavy-baryon ChPT (HBChPT) 
[g], 0], in which an expansion in inverse powers of the nucleon mass is performed in the 
Lagrangian. Later, several manifestly Lorentz-invariant renormalization schemes have been 
developed that also result in a proper power counting jl, 0, [l^, 11, 12, 13], with the infrared 
(IR) regularization of Ref. jol the most commonly used scheme. In its original formulation 



IR regularization is applied to one-loop intgrals, while the reformulated version of Ref. [14 



is also applicable to multi-loop diagrams [l5|. A different generalization of IR regularization 



to two-loop diagrams was suggested in Ref. [16 



Calculations at the two-loop level can be used to test the convergence behavior of the 
chiral expansion, which, compared to the purely mesonic part of the theory, seems to be 
slower in the baryonic sector. This is also of interest to lattice calculations, where chiral ex- 
trapolations are performed to obtain quantities at physical pion masses (see, e.g., Ref. 17]). 
The nucleon mass at order 0{q^) has been analyzed in the framework of HBChPT [l8|, in- 
cluding the evaluation of two-loop diagrams. In Ref. [l^ renormalization group techniques 
were used to determine the leading nonanalytic contributions to the nucleon axial-vector 
coupling constant qa at the two-loop level, which are independent of the applied renormal- 
ization scheme. To the best of our knowledge, however, no complete two-loop calculation in 
a manifestly Lorentz-invariant formulation of baryon ChPT (BChPT) has been performed 
so far. The determination of the nucleon mass up to a given order is one of the simplest 
calculations that can be performed in BChPT up to that order. This makes it the ideal 
physical quantity to perform a complete and consistent calculation at the two-loop level. 
In this Letter we present the results of such a calculation up to and including order 0{q^) 
using the reformulated infrared renormalization. 



II. LAGRANGIAN AND POWER COUNTING 

The effective Lagrangian relevant for the calculation of the nucleon mass to order 0{q^) 
is given by the sum of a purely mesonic and a one-nucleon part, 

^cS = £2 + A + ^nN + ^wN + ^wN + ^^wN + ^nN + ^wN + ' ' ' ■ (l) 

The purely mesonic Lagrangian at order 0{q^) is given in Ref. [2|. Reference contains 
the mesonic Lagrangian at order 0{q'^) as well as the lowest-order and next-to- leading-order 
nucleonic Lagrangians. The Lagrangians of the nucleon sector at order 0{q^) and 0{q'^) can 
be found in Refs. 2^, 21|. The complete Lagrangians at order 0(g^)and 0{q^) have not yet 



been constructed. However, up to the order we are considering, the nucleon mass does not 
receive any contributions from the Lagrangian at order O(g^), since only even powers of the 
pion mass can be generated from contact terms. The contributions from the Lagrangian at 
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FIG. 1: One-loop diagrams contributing to the nucleon mass up to order 0{q^). 



order 0{q^) will be of the form giM^, where gi denotes a linear combination of low-energy 
coupling constants (LECs) from C 



(6) 
■kN- 



We use the following standard power counting 2^, 23 1: Each loop integration in n dimen- 
sions is counted as q^, a pion propagator as a nucleon propagator as and vertices 
derived from £j and C^^]^ as and q\ respectively. 



III. RESULTS AND DISCUSSION 

Figure [1] shows the one-loop diagrams which generate nonvanishing contributions to the 
nucleon mass up to order 0{q^) . Diagrams with mass insertions in the nucleon propagator are 
taken into account by shifting the mass in the undressed nucleon propagator [9|. Diagrams 
(a) and (d) are of order 0{q^) and 0{q^), respectively. Diagrams (b) and (c) are of order 
0{q^), while diagram (e) counts as 0{q^). The relevant two-loop diagrams are shown in 
Fig. O According to the power counting there are additional diagrams at the given order, 
which however give vanishing contributions to the nucleon mass up to this order. An example 
would be diagram [2] (c) with one first-order vertex replaced by a second-order one. 

The renormalization of one- and two-loop integrals is performed using the reformulated 



infrared renormalization of Ref . [ij] . Details of the renormalization procedure will be given 
in a forthcoming publication j23 |. 

Since the perturbative expansion is performed around a ground state realized in the 
Nambu-Goldstone mode, the quark mass expansion of physical quantities contains analytic 
as well as nonanalytic terms |25|. The chiral expansion of the nucleon mass up to order 
0{q^) reads 

rriN = m + fciM^ + k2M^ + k^M'^ In — + k^M"^ + k^M^ In — + keM^ 

/i n 

+k7M^ In^ — -I- ksM^ In — + kgM^ , (2) 

where = 2Brh is the leading-order expression for the pion mass squared in terms of 
the average light quark mass, m = (m^ + md)/2 [2^, and m stands for the nucleon mass in 
the chiral limit. Combining the contributions from contact interactions with the one- and 
two-loop results the coefficients ki are given by 



(a) 



(b) 



(c) 



(d) 



-®—®—(^ — + — ®—(^)—(iy 



(e) 



(f) (g) 



(h) 



(i) 



-<1^ (i>- 



-<D ^ 



(j) 



(k) 



(l>- 



(1) 



FIG. 2: Two-loop diagrams contributing to the nucleon mass up to order 0{q^ 
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+ [3 - 11527r2F2ciC2 + 11527r2F'm 63 + 3207r2F2m eigl . (3) 

To simplify the notation we use 

ei = 16638 + 2eii5 + Scne , 

62 = 2ei4 + 2ei9 - 635 - 4e38 , 

63 = ei5 + 620 + 635 

for combinations of fourth-order baryonic LECs, while gi denotes a combination of LECs 
from the Lagrangian at order 0{q^). 

In general, the expressions of the coefficients in the chiral expansion of a physical quantity 
differ in various renormalization schemes, since analytic contributions can be absorbed by 
redefining LECs. However, this is not possible for the leading nonanalytic terms, which 
therefore have to agree in all renormalization schemes. Comparing our result with the 



HBChPT calculation of [18|, we see that the expressions for the coefficients ^2, ^3, and 
/cs agree as expected. At order 0{q^) also the coefficient kj has to be the same in all 
renormalization schemes. Note that, while k^M^ and fcgM^ln^ are nonanalytic in the 
quark masses, the algebraic form of the coefficients kg and kg are renormalization scheme 
dependent. This is due to the different treatment of one-loop diagrams. The counterterms 
for one-loop subdiagrams depend on the renormalization scheme and produce nonanalytic 
terms proportional to and In ^ when used as vertices in counterterm diagrams. 

The numerical contributions from higher-order terms cannot be determined so far, since 
most expressions in Eq. ([3]) contain unknown LECs from the Lagrangians of order 0{q'^) 
and higher. Only the coefficient /C5 is free of these higher-order LECs and is given in terms 
of the axial-vector coupling constant qa and F. While the values for both qa and F should 
be taken in the chiral limit, we evaluate k^ using the physical values qa = 1.2695(29) [2^ 
and = 92.42(26) MeV in order to get an estimate of higher-order corrections. Setting 
/i = niN, rriN = {rrip + m„)/2 = 938.92 MeV, and M = = 139.57 MeV we obtain 

k^M^ ln(M/mAr) = —3.8 MeV. This amounts to approximately 25% of the leading nonana- 
lytic contribution at one-loop order, k2M^. For a discussion of the importance of the terms 
at order 0{q^) at unphysical quark masses as used in lattice extrapolations see Ref. 271 . 



IV. SUMMARY 

Using the reformulated infrared regularization [1^] we have calculated the nucleon mass 
up to and including order 0{q^). This is the first complete two-loop calculation in manifestly 
Lorentz-invariant baryon chiral perturbation theory. The applied renormalization scheme 
preserves the standard power counting and respects all symmetries. Our results for the 
renormalization scheme independent terms agree with the HBChPT results of Ref. [1^] . 
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